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I. INTRODUCTION 



The coherent states are a fundamental tool in quantum mechanics. They were initially 
envisioned by Schrodinger^ as minimum uncertainty Gaussian states whose dynamics has 
maximum similarity to the classical oscillator. The interest in these states, which are related 
to the Weyl-Heisenberg algebra, was again aroused by the work of Glauber-'"^ Klauder''"^ 
and Sudarshan'', with the emergence of the first applications directed to quantum optics and 
semiclassical approximations. 

The generalization of the coherent states to arbitrary Lie groups was initially studied 
by Klauder, but a complete and detailed definition of the generalized coherent states was 
only later developed by Perelomov' and Gilmore'". In this way, important properties of 
the harmonic-oscillator coherent states were extended to other algebras of physical interest, 
resulting in an ideal setting for the study of the quantum-classical correspondence^^. 

Although the semiclassical methods also have their origin in the early days of quantum 
mechanics, the first calculations concerning semiclassical propagators in the representation 
of harmonic-oscillator coherent states were only developed in the last decades of the past 
century^-"^ ', as an attempt to transpose the practical difficulties of semiclassical propagators 
in the coordinate and momentum representations. 

However, similarly to the very definitions of coherent states, the generalization of semi- 
classical propagators for arbitrary groups occurred only in the decades following the semi- 
classical approximations with Gaussian states. In particular, aiming at the treatment of 
systems with spin, various correct derivations of the SU(2) semiclassical propagator ^''"^^ 
were obtained independently. Amidst these different methods, Kochetov^'*'-" generalized his 
results to arbitrary Lie groups. 

The study of semiclassical propagators is motivated mainly by the computational gain in 
relation to exact quantum calculations, considering practical applications in systems with 
many degrees of freedom or large quantum numbers, in which cases it is believed that the 
semiclassical approximation is efficient and accurate. From the less pragmatic point of view, 
the use of classical quantities for the calculation of quantum corrections allows us a greater 
understanding of the quantum-classical correspondence principle, since it establishes more 
precisely the boundaries and intersections between the two theories. 

We propose in this paper an independent approach for the derivation of the semiclassical 
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propagator in arbitrary coherent state representations. However, to make this text clearer 
and more accessible, here we exemplify our method only for the coherent states related to the 
fully symmetric representations of the SU(n) group^^ These coherent states are appropriate 
for the treatment of bosonic systems with a finite number of modes and conservation of 
the total particle number. Therefore, the SU(n) coherent states are important for several 
physical applications, as for example, the research on Bose-Einstein condensates in multi-well 
traps-^^"^*^. 

The major technical difficulty encountered in the derivation of semiclassical propagators 
is related to the explicit calculation of path integrals for quadratic approximations of the 
action functional. The topology and dimensionality of the classical phase space can make 
this task quite complex. Recently, Braun and Garg-'^''^*^ developed a simple and transparent 
solution to the semiclassical propagator for an arbitrary number of degrees of freedom, but 
considering only the direct product of many harmonic-oscillator or SU(2) coherent states. 
Using a slightly modified form of their previous findings and a suitable transformation of 
variables, we extend the path integral calculation to the SU(r;,) coherent states. 

The remainder of this paper is divided as follows: in large part of section II we reproduce 
the main results of Braun and Garg"' ' concerning the semiclassical propagator for multidi- 
mensional harmonic-oscillator coherent states. At the end of this section, we reformulate 
the identity that relates the path integral resulting from the semiclassical approximation to 
dynamical quantities calculated on the classical trajectory. This fundamental result is used 
in section III, where we make a detailed derivation of the semiclassical propagator in the 
SU(n) coherent state representation. Finally, the section IV is intended to our concluding 
remarks. 



II. SEMICLASSICAL PROPAGATOR FOR (n - 1) HARMONIC MODES 
A. Harmonic-oscillator coherent states 

The coherent state for [n — 1) harmonic modes is given by: 

\z)=e~^ Y] TT r — ^ |mi,m2,...,m„„i); (1) 

mi,m2,...,m„_i=0 \j=l V J / 

where {|mi,m2, . . . ,m„_i)} is the usual basis of the bosonic Fock space B*^"^ for (n — 1) 
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modes, such that rrij is the occupation in the j-th mode. The complex vector z = 
{zi, Z2, . . . , Zn-i)'^ , with {n — 1) entries, parameterizes the entire set of coherent states. 
Also, according to the adopted notation, the juxtaposition of two vectors a and b represents 
the matrix product ab = aibi + 02^2 + • • • + an_ib„_i. 

The coherent states (1) constitute an overcomplete set in Hilbert space, which enables us 
to write the following diagonal resolution of the identity in B"^^: 

df,iz*,z)\z){z\=l; dfj.iz*,z) = ll^; (2) 

1=1 ^ 

where d'^Zj = dqjdpj = — with qj = Re (zj) and pj = Im (zj). 



B. Coherent state propagator 

The quantum propagator in the coherent state representation is defined as the matrix 
element of the time evolution operator between the initial coherent state l^j) and the final 
coherent state \zf), assuming a propagation period r: 

K{z},zf,T) = {zf\e-'''^\zi). (3) 

In the above definition, for simplicity, we assumed h = 1 and the temporal independence 
of the Hamiltonian H. However, the main results presented in this work, including the semi- 
classical approximation of the propagator, are also valid for Hamiltonians with explicit time 
dependence"' . Expanding the propagator (3) to second order around the classical trajectory, 
which is defined later, we obtain the following form for its semiclassical approximation: 

Ks,{z}, z,; r) = e^^^-^d^'I'+l^/l'^K.e^l^, z,; r); (4) 

where Sc is the action functional S calculated on the classical trajectory, and Kred is the 
reduced propagator, which introduces quantum corrections. The classical trajectories, by 
definition, are the extremes of the action functional: 

r 

iSiz},Zi;T) = J !^^i-2z-zz)-t'Hiz,z)^dt+^[z}ziT) + ziO)z^]. (5) 
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In the previous equation, to emphasize the independence of the variables z and z* ^ we 
made the change of notation z* — z. Also in (5), we defined the function 'Hl'ZjZ) = 
^\J*\zf ' which represents the effective classical Hamiltonian. The extremization of the action 
functional implies the following classical equations of motion'^'^: 

" = -^^' ' = '^- 

According to (3), the correct classical trajectory must respect the boundary conditions 
z{0) = Zi and z(r) = z^ Generally, if we consider the variables 1 as the complex conjugates 
of z, then the boundary conditions would make the system of differential equations (6) 
overdetermined, since the two vector equations would become redundant. Therefore, for 
arbitrary classical trajectories, is not the complex conjugate of z{t), henceforth denoted 
by z*{t), except for a specific boundary condition in which accidentally I'(O) = z*. 

The reduced propagator corresponds to an infinite dimensional integral of the second 
variation of the action 



Kred{z}, Zi] t) = D^ij], rj) exp I > ; Dfi{r], ij) = lirn^ JJ 



n-l j2 k' 
TT 



(7) 



The integration variables i]^ = z^ — z\ and if =l}^ — 'z'l represent the departure from 
the classical trajectory, denoted by z^ and zj, at the time tk = ke, for k = 0, 1, 2, . . . , M 
and e = jj. These variables are subject to the boundary conditions 77° = r]{0) = and 
_ fj(^f"j = 0^ also in agreement with the propagator (3). 

By linearization of the equations (6), we can also obtain the classical equations of motion 
for deviations from the classical trajectory. In this way we define the matrices Rj^, which 
are essential for later results^'^ 




azaz dz 





(8) 



dz^ dzdz 

The quantities in the previous equation are calculated on the classical trajectory. Notice 
that 6z and d'z, which represent the dynamical variables of equation (8), deserve different 
notation from the integration variables r]^ and tJ^ due to their fundamentally unrelated 
purposes and definitions. 
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After a lengthy derivation, the result of the Gaussian integral (7) can be written in terms 
of dynamical quantities calculated on the classical trajectory^*^: 



Kred{z}, Zi] T) = e 



f Tr[R22it)]dt 



'det 



^(0) 



(9) 



The matrix represents the response of the "free" trajectory extreme z{0) to a change 



in the boundary condition z{t) = zj. Employing the functional (5), we can show the identity 
relating the second factor of (9) to a second derivative of the action. 



9z(t) dzidzj 



C. Continuous limit of 

In the original text of Braun and Garg'^'', only the time discretized form of 6^Sc was 
employed in order to obtain the result (9). In this approach, the Gaussian integrals in 
(7) are calculated sequentially, before the application of the continuous-time limit {e — 0, 
M — )■ oo and eM = r). However, we can write the continuous form of directly from 
(5): 



r 

ib'^Sc = J dt[(j]r] - fjrj) - r]R2ir] - + rjRiiV + vRi2v} ■ 



Substituting (10) in (7), we obtain a new expression for the reduced propagator: 



(10) 



Kred{z}, Zi] ^) = y DiJ.{r], rj) exp | ^ j dt [(7777 - rjrj) - r]R2iri - 2r]R22'n + vRuv] | ; (H) 

where we used the identity R22 = —Rii- Comparing the representations (9) and (11) of 
Kred, we can write a formal identity between the path integral of d'^Sc in its continuous form 
and the corresponding solution described by classical quantities: 



^ f dt[{7jr]-riri)-r]R2iV-27jR22V+vRi2rj\ ^ f TT[R22(t)]dt 

D^{rj,ri) e « =e « Wdet 



&z{Q) 
Mr) 



(12) 



We must always remember that the previous expression is associated with the boundary 
conditions 77(0) = and fji^T) = 0. The equation (12) is the main result of this section, since 
its generalization to SU(n) coherent states is almost straightforward, as shown below. 
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III. PATH INTEGRAL FOR SU(n) COHERENT STATES 



E (^rw^) in»r)^=^. (13) 

mi+m2+...+mn=N 



A. SU(n) coherent states 

The coherent states related to the fully symmetric irreducible representations of SU(n) 
are given by^-^: 

N\ V (t]w"''^ |mi,m2, ■ ■ .,m n) 
mi\m2\...mn\) \^}^ ^ ) {l+w*wy 

Here again {|mi, 777,2, • • • , fn-n)} is the usual basis of the bosonic Fock space for n modes 
and particles, such that r?7j is the occupation of the j-th mode. Similarly to the vector z 
defined in equation (1), the number of entries of the complex vector w = {wi, W2, ■ ■ ■ , Wn-i) 
is (77 — 1). However, the number of bosonic modes is now n, but with occupations restricted 

n 

by ^ nij = N. Note that the total number of particles is the only index needed to 
uniquely specify a fully symmetric irreducible representation of SU(?7). 

The SU(n) coherent states are suitable for the treatment of bosonic systems with a finite 
number of modes and conservation of the total particle number. Therefore, henceforth we 
consider a system with a fixed number of bosons, so that the evolution of its state is restricted 
to only one SU(?7) irreducible representation. Thus, we also simplify the notation by making 
|A^; w) \w). 

Just as in the harmonic-oscillator case, the coherent states in (13) are normalized, but 
they are not orthogonal. The overlap between two SU(?7) coherent states is given by: 

{w'\w) = ^ (14) 

(1 + WW) 2 (1 + W*W) 2 



Using the SU(?7) coherent states, we can write a diagonal resolution for the identity in 



I df,{w\w)\w){w\ = l; df^{w\w) = ^^^^^^lld'w,. (15) 

Here again we use the notation d'^wj = dxjdyj = 'tHj^Hj.^ where Xj = Re (7^^) and 
yj = lm{wj). Notice that we can deform the integration domain by making wj and w* 
independent variables, which amounts to the complexification of xj and yj. We use this 



procedure in the following subsections to replace the complex plane by a surface of integra- 
tion in for each value of j. Also note that the normalization constant in (15) can be 
divided into the factors a{n) = ^""ij' , which is independent of the total particle number, 
and dim(B^) = ^^^^Iiy , the dimension of accessible Hilbert space. 

The choice of parameterization defined in (13) for the SU(n) coherent states is quite 
appropriate, since the integration domain in (15) is the same of equation (2). Moreover, 
except for the normalization factor, the coherent states are entire functions of the variables 
w. 

B. SU(n) coherent state propagator 

Similarly to equation (3), we define the propagator in the SU(n) coherent state repre- 
sentation as the transition probability between the initial coherent state \wi) and the final 
coherent state \wf) after a time interval r: 

K{w},Wi;T) = (wfle-'^'^lwi). (16) 

For simplicity, we have assumed again a time- independent Hamiltonian H. To write down 
a path integral representation for K{w*j-,Wi]T), we factorize the evolution operator into M 
identical propagation subintervals of length e = j^. Then, we introduce an identity of the 
form (15) between each pair of factors in the discretized evolution operator: 

M 

JJ(mJ^>-*^"|m;^-^). (17) 

i=i 

In the previous equation we define the integration variables w'^ = w{tk) and = w{tk), 
with tfc = ke for k = 1,2, ... , (M — 1). We also extend this notation to the propagator 
boundary conditions by setting = w{0) = Wi and w'^ = w{t) = w*j. Note that we 
employ the notation w*{t) — )■ w{t) since the beginning of this section, already evidencing 
the future duplication of the classical phase space. However, and need to be considered 
as independent integration variables only after the semiclassical approximation, when the 
integration domain is deformed to a surface in the doubled phase space which contains the 
k-i\i point of the discretized classical trajectory. Nonetheless, to maintain the consistency of 
notation in this procedure, notice that the identities nP = w^* = w* and w'^'^ = w^^* = Wf 

8 



K{w},Wi; t) 



M-l 



must always be preserved {w^ and can not be made independent of and w^), since 
these quantities do not stand for integration variables, which can be designated to follow 
the classical trajectory in the doubled phase space. 

Considering that at the some point of the calculations we will take the limit M — j- oo, 
we can expand each factor in the integrand of (17) to first order in e: 



(18) 



where "H 



{w^*\H\wi-^) 



= ^-^^jTj— ny^ , for j = 1,2,..., M. Then, substituting the expression (18) in 



(17), we can readily obtain: 



K{w*j-,Wi;T) 



k=l 



exp 



■ M 

El- 

.i=l 



M 



(19) 



Employing the equation (14), the argument of the exponential in (19) takes the following 
form: 



iSd 



M 

E 

i=i 
i 



-Ln 



1 + w^w 



_(1 +Ww^){l +W^^wi'^) 



(20) 



5d-yLn [{l + \w}\^){l 



\Wi 



In order to achieve the correct continuous limit of Sd-, it is imperative to note that generally 



w{t) 7^ w = Wf and w{0) ^ 



w*, thus ensuring the continuity of the trajectories in 



the independent variables w{t) and w{t). Although explicitly depends on w^^ and uf, as 
we see in the first line of (20), it is possible to analytically isolate all dependence on these 
quantities in a single term, as shown in the last line of (20). The term Sd represents the 
discretized form of the action functional '' , whose extremization correctly gives the classical 
equations of motion in the continuous-time limit. 

Also in this limit, we can assume that consecutive values of w{tfc) and w{tk) are very close, 
forming a continuous trajectory. Therefore, according to (14), we can write the following 
identity, which is valid up to first order in \w{tk) — 'w{tk)\ and \w{tk) — w{tk)\: 



Ln{w*{tk)\w{tk-i)) 



N [w{tk) - w{tk-i)]w{tk-i) - w{tk)[w{tk) - w(tfc-i)] 
2 1 + W{tk)w{tk-i) 



(21) 
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Adding and subtracting the appropriate terms dependent on w{t) and w{0), we can easily 
take the continuous hmit of the action functional Sd with aid of expression (21): 



r 

tSiw},w,;T)= [dt\^ WW -WW _ }N^^^ ^ _ ^ ^^2) 

■' J [21 + WW } 2 ^ ■' 



In the previous equation we also defined the effective Hamiltonian T-Linj^w) = , 
which represents the continuous limit of 'Hjj-i. The extremization of S, for fixed boundary 
conditions and propagation period, provides the classical equations of motion. So, by making 
6S = 0, we obtain'^-^: 

w = - — (1 + WW) [l+w xw) — = -ir.—; 

iV OW OW (nn\ 

^ ^ _ W. - N --^^ 

W = — I 1 + WW) [1 + W X W) — — = l^T^- 

N OW OW 

Here we should remember again that the classical trajectory is contained in a doubled 

phase space, since w{t) and w{t) must be treated as independent variables, preventing the 

overdetermination of the equations of motion caused by the boundary conditions w{0) = Wi 

and w[t) = Wj. 

Analogously to (8), we can write equations of motion for deviations from the classical 
trajectory. Thus, we redefine the matrices Rjk for our new dynamical variables: 




(24) 



\~dH] -iA. l^au] \ Sw 

dw I dw J dw I 9w J I I 

dw I dw J dw I dw j / \ 

Notice that the identity R22 = —Ru is no longer valid. 

C. Semiclassical approximation 

Substituting the definition (20) in equation (19) and then taking the continuous limit, 
we obtain a new expression for the SU(n) propagator: 



K{w},Wi;T) = j Dfi{w,w) exp i^iS{w},Wi;T) - ^Ln [(1 + |«;/|^)(1 



w.n\ . (25) 
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The quantity Djj,{w, w) symbolizes the following infinite product of measures in the dou- 
bled phase space: 



Dii(w,w) = lim If 



M-1 

A/-i>oo 



.(.)din.(B~)||^, - 



(26) 



fc=i 

Now, we proceed to the semiclassical approximation by expanding the action up to second 
order around the classical trajectory, which we denote by Wcif) and Wc{t): 

S^S, + ^6^S,; 6S, = 0. (27) 

Substituting the above approximation in the propagator (25), we obtain the semiclassical 
propagator in the SU(n) coherent state representation: 



Therefore, it remains only to calculate the reduced propagator, defined by: 



(28) 



Krediw}, w,; t) = J Dfi{r], rj) e^'"^^- (29) 

where we have introduced new integration variables rj{t) = w{t) — Wc{t) and f]{t) = w{t) — 
Wc{t), which represent deviations from the classical trajectory. Note that these new variables 
are subject to the boundary conditions r]{0) = r/(r) = 0. Expanding the action (22), we get 
the following quadratic terms in r] and f] after some integrations by parts: 



T 

id'^Sc = J dt {riei]-r}Qr] + r]Ar] + 2r]Bf] + r}Cf]} . (30) 



The matrices and O are the inverses of S and H, respectively. The explicit form of 
these matrices is given by: 



e = ^ i^ + ^^n-^^^ = e^. (31) 

The other quantities introduced in equation (30) are defined as follows: 
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^ 2{Ww)w X w — {1 +Ww){w X w + w X w) i d^H 

(1 + wwy N dw'^ ' 

^ 1 {Ww — Ww)[2w X w — {1 + Ww)l] + (1 + ww)(W x w —W x w) i d^H _ ^^^-^ 

2 (l + ww)^ N dwdw^ 

^ (1 + ww){w X w + w X w) — 2(Ww)w X w i d^T-L 

(1 + Ww)'^ N dw"^ ' 

Now we make a further change of integration variables, in order to compare the path 
integral (29) with the canonical result (12). To this end, we define the following matrices: 



„ /—IWXw + ^/T+ww[{ww)l—wxw] —T , . 

Q = \ N — — — : r = Q ; (33) 

ww[l + ww) 

such that = Q and Q =0. Then, applying the linear transformations v = Qrj and 
V = Qrj in (30), we get: 

r 

iS'^Sc = J dt I (z/77 - z7z>) + uAu + 2iyBV + VCV^ . (34) 







Note that, in this new set of variables, has the same form previously found in (10). 
The matrices above are related to variables of interest w and w as follows: 



B = Q 'bQ '-1(qQ 'g) ; (35) 



A = Q 'AQ-'; 

~ 2 

C = Q-^CQ'\ 

In the previous definitions we used the inverses of the matrices Q and Q, which can be 
readily found: 

^ 1 1 + WW { \{ww)\ — W X w]^ f—-i\T 

= \wxw+ ^^ fY^ ^^ = [Q ) • (36) 

y/Nww I V 1 + WW J ^ ^ 

Under the same transformation applied to (34), the path measure (26) takes the following 
form: 



Daily, u) = lim TT 

M-^no J- -I- 



M-i-oo 

k=l 



(N + n-l)\''T^d^iy: 
— 11 



(37) 



Therefore, Dfi^u, v) takes the same form found in (7) only when N ^ n. That is, the 
transformation given by (33) produces the desired result only for large number of particles 
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(macroscopic limit). However, this is precisely the expected situation for the effectiveness 
of the semiclassical approximation, considering the computational gain and the quantitative 
agreement with respect to exact quantum results. 

In order to go from equation (26) to (37) we employed the explicit expression for det Q, 
which is trivially related to det G. The determinant of can be easily calculated from the 
identity: 

det [(1 + ww)l - wxw] = {1 + wwy-'^] for n>2. (38) 

Substituting (34) and (37) in (29), under the additional condition N ^ n, we can recast 
the reduced propagator in a similar way to the equation (11): 



red\ 



{w*f,Wi]T) = / Dfi{iy,i^) exp 



dt 



[vv — uu) + uAu + 2vBv + vCv 



(39) 



The variables v{t) and z7(f) can be interpreted as deviations from the classical trajectory 
described in a new set of dynamical variables, such that z/(t) = v{t) — Vc{t) and V{t) = 
v{t) — Vc(t). With this interpretation in mind, we can directly compare the equations (12) 
and (39), allowing us to solve the integral of the reduced propagator in an immediate way: 



(40) 



Kred{w*f,Wi]T) = exp < Tt B{t) 




However, the above result is not yet the desired solution, because we still need to return 
to the original variables w and w. Comparing the transformation between the variables r] (jj) 
and u (y) with their respective definitions, it becomes evident the relation v{t) = Q{t)w{t) 
(y{t) = Q{t)w(t)). Therefore, we can write the following identity: 



det 



dv{0) 
dv{T) 



det 



ow[t) 



1 + w{t)w{t) 



det 



dwir) 



(41) 



In the last equation we used the expression for the determinant of Q{t), which can be 
easily calculated from (38). Furthermore, considering the quantities defined in (24), (31), 
(32) and (35), we can also obtain the relation: 



Tr 



B{t) 



Tr 



e {t)B{t) =-1i[Ru{t)-R22{t)]. 



(42) 
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At last, substituting (41) and (42) into (40), we find the final expression for the semiclas- 
sical propagator of bosons in n modes: 



1 + w{t)w{t) 
_l + w{0)w{0)_ 



det 



dw{0) 
dw{T) 



(43) 



In the previous equation, we introduced the correction term to the classical action in this 
quantization scheme: 



7 = i/lV 





dw \ dw J dw \ dw J 



dt. 



(44) 



In comparison with the results for the harmonic-oscillator coherent states, we can also 
relate the factor under the square root symbol in (43) with a second derivative of the classical 
action: 



1 + w{t)w{t) 
l + w{0)w{0) 



det 



dw{0) 
dwir) 



[1 + WfW{T)]'^[l + w{0)wi]'^ det 



(45) 



IV. CONCLUSION 



We presented in detail a derivation for the semiclassical propagator in the SU(n) coherent 
state representation based on previous findings for the multidimensional harmonic-oscillator 
coherent states. The main advantage of our method is that it can easily be extended to 
other groups of physical interest. The transformation given by (33), which allows a direct 
comparison between the SU(r2) path integral and the results for the harmonic-oscillator 
coherent states, enables the immediate calculation of the semiclassical propagator in a phase 
space with more complex topology. This simple procedure can be readily generalized to a 
wide range of dynamical groups, thus avoiding an additional laborious derivation of the 
semiclassical propagator path integral for each new choice of coherent state representation, 
as usually done in previous publications-^''^\ 

The equation (43) has direct application to the dynamical analysis of bosonic particles 
in n modes. It is expected that the SU(?7,) semiclassical propagator provides accurate results, 
when compared to exact quantum calculations, for sufficiently large values of A^. However, 
we can safely suppose that the approximation prescribed by AT^c exceeds significantly the 
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accuracy of mean-field methods, which are usually employed in the treatment of bosonic 
dynamics, even for a relatively small number of particles'^^ . We also expect the semi classical 
propagation to be computationally more efficient, in relation to similar quantum calculations, 
for systems with many degrees of freedom (large values of n). In this case, the computational 
cost of calculation for the involved classical quantities grows slower with the number of modes 
than their quantum counterparts. 

ACKNOWLEDGMENTS 

We acknowledge the financial support from CNPq and FAPESP, under grants No. 
2008/09491-9 and 2009/11032-5. 

REFERENCES 

^E. Schrodinger, Naturwiss. 14, 664 (1926). 
2R. J. Glauber, Phys. Rev. Lett. 10, 84 (1963). 
3R. J. Glauber, Phys. Rev. 130, 2529 (1963). 
^R. J. Glauber, Phys. Rev. 131, 2766 (1963). 
5J. R. Klauder, Ann. Phys. (N.Y.) 11, 123 (1960). 
'^J. R. Klauder, J. Math. Phys. 4, 1055 (1963). 
7j. R. Klauder, J. Math. Phys. 4, 1058 (1963). 
^E. C. G. Sudarshan, Phys. Rev. Lett. 10, 277 (1963). 
^A. M. Perelomov, Commun. Math. Phys. 26, 222 (1972). 
ioR. Gilmore, Ann. Phys. (N.Y.) 74, 391 (1972). 

i^W.-M. Zhang, D. H. Feng, and R. Gilmore, Rev. Mod. Phys. 62, 867 (1990). 

i^Y. Weissman, J. Chem. Phys. 76, 4067 (1982). 

i^Y. Weissman, J. Phys. A: Math. Gen. 16, 2693 (1983). 

i^D. Huber and E. J. Heller, J. Chem. Phys. 87, 5302 (1987). 

^^D. Huber, E. J. Heller, and R. G. Littlejohn, J. Chem. Phys. 89, 2003 (1988). 

^^R. G. Solari, J. Math. Phys. 28, 1097 (1987). 

i^V. R. Vieira and P. D. Sacramento, Nucl. Phys. B 448, 331 (1995). 
i^E. A. Kochetov, J. Math. Phys. 36, 4667 (1995). 

15 



i^E. A. Kochetov, J. Math. Phys. 36, 1666 (1995). 

20E. A. Kochetov, J. Phys. A: Math. Gen. 31, 4473 (1998). 

2iR. Gilmore, C. M. Bowden, and L. M. Narducci, Phys. Rev. A 12, 1019 (1975). 

22T. F. Viscondi, K. Furuya, and M. C. de Ohveira, Ann. Phys. (N.Y.) 324, 1837 (2009). 

23F. Trimborn, D. Witthaut, and H. J. Korsch, Phys. Rev. A 77, 043631 (2008). 

24F. Trimborn, D. Witthaut, and H. J. Korsch, Phys. Rev. A 79, 013608 (2009). 

25c. Braun and A. Garg, J. Math. Phys. 48, 032104 (2007). 

26c. Braun and A. Garg, J. Math. Phys. 48, 102104 (2007). 

2'^M. Baranger, M. A. M. de Aguiar, F. Keck, H. J. Korsch, and B. Schellhaafi, J. Phys. A: 

Math. Gen. 34, 7227 (2001). 
2^ According to the adopted notation, considering a scalar function /(a) and a vector a, we 



have that represents the vector whose elements are given by 



dfja) 
da 



^ dm 

dai ' 



for 



j = l,2,...,(n-l). 

2^Following the adopted notation, considering a scalar function /(a, b) and two arbitrary 
vectors a and b, we have that represents the matrix whose elements are given by 
9^] - for J, A; = 1, 2, . . . , (n - 1). 



jj. dajdb^.- 



dadb 

■^^For two vector quantities a and b, we denote by |^ the matrix whose elements follow from 

da^ 

^^The subindex d indicates the time discretization of dynamical quantities. Therefore, the 



[S],. = fe'fo^^'^ = l'2,...,(n-1). 



omission of this index points out the continuous limit of the corresponding functionals. 
^^Here we introduce the notation for the dyadic product between vectors. Considering two 

arbitrary vectors a and b, the matrix elements of a x 6 are given by (a x b)jk = djbk- 
33A. D. Ribeiro, M. A. M. de Aguiar, and A. F. R. de Toledo Piza, J. Phys. A: Math. Gen. 

39, 3085 (2006). 
■^^T. F. Viscondi and M. A. M. de Aguiar, in preparation. 



16 



